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Coefficient matrix (3 x 3) 
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Symbol 

p,q>r 


V, v : 


x,y,z 


Definition 

Gun angular rates about the 
aircraft x, y, z axes 

Unit vector in the direction 
of the target 

Bullet time of flight 
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Target velocity 

Target acceleration 
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I. INTRODUCTION 

With the installation of digital computers on fighter 
aircraft comes a computational capability that allows 
greater accuracy in the solution of fire control problems 
as well as greater versatility in implementing the solutions. 

This report is concerned with the feasibility of imple¬ 
menting the solution of the l_ead computing optical sight 
(LCOS) in an on-board digital computer. For purposes of 
this report this sight will be called Air Force Academy LCOS 
(AFALCOS). It is shown that the lead angle equation required 
for reticle displacement on the heads-up display (HUD) can 
readily be solved on the digital computer. The solution of 
the lead angle equation is compared to the steady state solu¬ 
tions from models of a present day advanced fighter LCOS, the 
Honeywell, Inc. digital LCOS, the Avionics Laboratory digital 
LCOS documented in reference [2J 1 , and Honeywell's historical 
tracer sight "HOTLINE." The AFALCOS is shown to be generally 
more accurate that the above LCOS's and in the steady state 
to compare very favorably with HOTLINE. 

It should be noted that comparing two generically differ¬ 
ent types of gunsights requires not only a steady state com¬ 
parison for accuracy but also a dynamic comparison in which 
the entire pilot, target, sight display, flight control system 



l Air-to-Air Gun Fire Control Equations for Digital 
Lead Computing Optical Sights, AFAL-TM-74-9-NVE-1, 
Apr 1974, Air Force Avionics Laboratory, Wright- 
Patterson AFB, Ohio. 
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and aerodynamic response of the aircraft are included. This 
complex interaction is highly coupled and the entire system 
needs to be tested before suitability of any one sight can 
be fully determined. In this report only the steady state 
comparison will be made. 

In Section II the problem definition is presented. The 
assumptions are stated, and the development of the equations 
is made in Section III. A discussion of the problem of 
implementation is made in Section IV. Numerical results 
demonstrating the feasibility of the digital LCOS are pre¬ 
sented in Section V as well as an accuracy comparison with 
the models stated above. Conclusions are stated in Section 
VI. 
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II. PROBLEM DEFINITION 


Definition of LCOS 

In the lead-computing gunsight, a prediction angle 
is continually computed and displayed on the Heads Up 
Display (HUD) as a reticle offset from a fixed reference 
on the gunsight. A typical reticle presentation on the 
HUD is shown in Figure 1. 


Zenith Tabs 


Center Pipper 
Roll Index 


Roll Tab 



Range Bar 


Figure 1. Typical Sight Reticle Display 


The prediction angle displaces the reticle in azimuth and 
elevation and is the angle between the present line of 
sight (if the pilot has the reticle superimposed on the 
target) and weapon line. The prediction angle accounts 
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for ballistic curvature, own notion* and predicted target 
motion. In the implementation of this system, several draw¬ 
backs are apparent. One is the necessity to predict future 
target motion so that the proper lead angle may be obtained. 
To make a prediction, some assumptions must be made concern¬ 
ing future target motion which are not always valid if the 
target is performing evasive maneuvers. Also, the reticle, 
if undamped, moves abruptly in the HUD due to abrupt control 
inputs which are common in combat. To overcome this move¬ 
ment, damping must be induced in the reticle, which has the 
adverse effect of making the reticle lag behind real time. 
Thus, time is needed for the on-board computer to settle to 
a solution. This settling may never occur if aircraft atti¬ 
tude and velocity are constantly changing. LCOS, therefore, 
has limitations which are undesirable. However, LCOS has 
proven itself worthy in close-in air combat and is a usable 
sight. 

The sight presently incorporates the use of a Lead Com¬ 
puting Gyroscope (LCG) which can be eliminated if the air¬ 
craft has either a strapdovm or ine-tial navigation system 
on board which will provide the aircraft body rate vector. 
The problem, then, lies in developing a digital program 
capaDle of taking the place of the present LCOS, including 
the lead computing gyroscope. This approach should yield a 
gunsight which is more accurate than the present one at a 
decreased initial cost and a reduced maintenance cost. 

10 . 
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III. DERIVATION OF LCOS EQUATIONS 
In this section the lead angle equations will be 
developed along with the time of flight equation. Results 
of reference [l] 1 are used extensively. 



Figure 2. Lead Angle Components 
From inspection of Figure 2, 

Ds ♦ Vyr f + (V T ♦ A g g z) T f Z /2 - D f c. 


For the small angles involved, 
(V + V ) c s v a + V F 


— ~ a — , m xr 

c "* V +V m 3 V +V b 
am am 


(Assumption). 


By the definitions of the four variables 

D f B < V a + V T f ‘ 

Substitute equations (2) and (3) into equation (1) 
to obtain 

l F-4fc Lead Computing Optical Sight System (LCOSS) 
Improvement Study, Report MDC A0610 1 June 1971, 

Me Donnell Aircraft Company, St. Louis, Missouri. 


( 1 ) 

( 2 ) 

(3) 
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s x (V T - V a a) ■ D (s x t) * D [t - (F . s)sj. 



The terra in brackets on the r.h.s. of Eq (10) was 

simplified to just F in previous LCOS developments, which 

led to an incorrect result in the effect of aircraft roll 

rate on the lead angle. 

• • 

NOTE: s' * F x s’ (derivative of a unit vector) 

cross both sides with s' 

• • • • • 
s' x s = s' x (F x s') and s x s = F • (8 . s) s 


and s x a = A + (F x a) = A + a. 


(ID 


where s x F = A 


(Assumption) 


»q (Assumption). 


D x a = aj 

Substitute equations (8), (10), and (11) into (7) to obtain 


'■ (v 7 ;)'' r ’ !'■ '•' 


\ a m / 


( 12 ) 


(s x V T ) T f /2 + Ag g (s x z) T f /2. 


* D [F - (8 . s)s] + V a^ V m~ V f^ (A + a r ) 


V ♦ v m 

a m 


+ (s x V T ) T f /2 + A r g (s x z) T f /2. 


(13) 


Collect terms in A: 


V f A = D8’ + J y V a a G + A^ g (s x z) T f /2 + (s x V T ) T-/2 
where 


Jv * ^m ~ ^f and 8’ = F - (8 . s) s. 


V ' + V m 
a m 


(14) 

( 15 ) 
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Target acceleration is the only component of lead which is 
not expressed directly in terms of quantities measurable on 
board the aircraft. The next objective is to relate target 
acceleration to those parameters which can be either effec¬ 
tively measured or estimated. 

From equation (9) , 

V T = V„ a + d (Ds) = V a + D? + Ds. (16) 

T a cfif a 

At this point, previous derivations of LCDS have assumed 
constant target and attacker speeds; however, this assump¬ 
tion is not necessary unless only the normal acceleration 


is measured. Proceeding without this acrumption, we can write 

V r = d (V a) + d (D i x s) + d (Ds). (17) 

T It a <Tt eft 

We now write Newton's second law for the attacker, 

J T+D + L + W = md (V_a). (18) 

dt a 

Rearranging (18) we have 

d (V o aj = (7 + E + D/m + g. (19) 

3t a 

The right hand side of (19) is arranged in the manner 
shown because an on-board accelerometer triad will 
measure only nongravitational accelerations. Let a^ 
denote the nongravitational accelerations, 
a N = (7 + D* + D/m. (20) 

Equation (19) becomes 

f t - 5 n * *• < 21 > 

Substituting Equation (21) into Equation (17) and taking 
the vector cross product of both sides of equation (17) 
results in 
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where F is the gunline unit vector. From Figure 3 it is 
clear that 

] 1c | = A and 

s*F«-lc = F + (F x X). 



Figure 3. Relationship of A, s, F, and X 

In the r.h.s. of equation (24) we can express s x a,, in 

Is 

terms of a^, F, and X, all of which can be measured or 
computed on board the attacker aircraft, 

s x a N - [F + (F x A)] x a„ 


= F * a N • (V*)b * (v S)I - 

Next we consider 3- Referring to Figure 4 we can write 


F s « - X 


(Assumption). 

s’ (Line of sight to target) 


F (Gunline) 


Figure 4. Geometrical Relationship Between 6, w, and A 
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We nov need to simplify the terms in equation (24) 
containing F which are 

T 

(D + T f D) [3 ■ (3 ♦ s)s] + 'j— D (3 • s)s x IT* 

• • 

The term [F - (F . s)?] can be written in matrix 
form as 


6’ * [F - (3 • s)s] 


’(1-sf) - s x s. 

_ ■ Si s 2 (1-sl) 


"Si S3 * S2 S3 


s 1 s 3 


* S 2 S3 


(1-sf) 


= [A] 3. 


The term s x 3 can be written in matrix form as 


*0 -s 3 

«T x F = S3 0 

- s 2 si 0 


0 -Si 


= [B] 3. 

Finally, we assume that 


F • s = p 


(Assumption) 


where p is the roll rate of the aircraft. 

This assumption is very reasonable since F is 

the angular rate vector of the line of sight and F • s' is 

simply the component of angular rate along the line of sight. 









The line of sight vector s’ is nearly aligned with the 
roll axis of the aircraft. 

The expression (29) can now be written 


(D + T f D)[A] 8 + y- D p [B] 8 

- (D + T f D){[A] + T f Dp [B]} i 

2(D+T f D) 

- (D+T f D) (A + bB} i. 

Substituting equation (33) and equation (28) into the lead 
angle equation Eq (24) results in 


V f X = (D+T f D) (A + bB) (w-X) + J y V a a Q * T _± s x a N 

2 


♦ (1-Ag)*f g s x z ♦ D ' S f [A] 8 • 

M 4 


The last two terms are neglected in current LCOS 
mechanization and will be neglected here since 
the terms (1-Ag) and |F| are generally very small. 

We desire a differential equation for the lead angle; 
therefore, we solve for X. Rearranging Eq (34) and 
neglecting the last two terms 


(A+bB}X = {A+bB)w - V f 


X + J v V a 


(D+T f D) 


(D+T f D) 


T 

+ A f 


s x a. 


2(D+T f D) 


(35) 







s&jgS&id^ 



To solve for A it is necessary to invert the matrix 


{A+bB}; however, this matrix is singular if the line of 

sight vector corresponds with the aircraft body axis. The 

♦ 

matrix {A+bB} multiplying \ in Eq (35) can be written 

{A+bB} * (1-s^) -( s ! s 2 + bs 3 ) C-s l s 3 + bs 2 ) 

(-SjS 2 +bs 3 ) (l-s|) - (s 2 s 3 + bSj) ( 3 ( 

-(sis 3 + bs 2 ) (-S2S3 + bs\) (l-s§) 

Note th*t if Si = l, s 2 = s 3 = 0, then the first row of 
(A+bB) is zero, and the matrix is singular. This problem 
can be circumvented by deleting the first component 
equation for X 3 since X 3 is not required for displaying 
the lead angle. However, before deleting the first row 
and column of {A+bB} it is necessary to first multiply the 
matrices {A+bB} u on the r.h.s. of Eq (35). This multi¬ 
plication introduced the effect of roll, p, into the. equations 
• • 

for X 2 and X 3 . The resulting equation is 


(1 si) 

(-s 2 s 3 + bs j) 


-(s 2 s 3 + bs 3 ) 
(1-si) 


* T t- 


2(D*T f D) 


_ _ V 

s x a N - _f_ 


D+T f D 


_ (-SjS 2 + bs 3 )ui + (l-s 2 )w 2 - (s 2 s 3 + bs 1 )oj 3 

2 

(“S 3 s 3 + bs 2 )wi + (- s 2 s 3 +bsi)a > 2 + (l“S 3 ")w 3 j 


& 


J V 
+ v a 


(D+T f D) 

























and divide both sides of equation (48) by T £ : 
R/l \ 2 - V /1 \ + fc.(V. + R)V_' 3 p /p 1- 0 


.(IV [V V a * « V p % P /%] 


1 - v__ + \/v_ 


4V LS CV a + R) . 




Letting 



(53) 

(54) 


Equation (54) is used to calculate time of flight T £ for 
the lead angle equation, equation (43). To obtain V £ , the 
average bullet speed with respect to the attacker, 



V £ = D/T f + R . 



V 

a 


(55) 


(56) 


Equation (56) is used to compute V £ for the lead angle 
equation, equation (43). 


Note in equations (53) and (54) that V ; hence, T- may be 
undefined if the argument of the square root is negative. 
This situation can occur for certain engagement conditions 
and must be checked for solution validity. 
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IV. IMPLEMENTATION OF LCOS EQUATIONS 


In practice, the z vector may not be available. The 

term containing z is very small in any case (ref 1); hence, 

•• 

this term is often neglected. The F vector can be obtained 

•• 

by differentiating F in equation (28); however, because 
differentiation is a noisy process and this term is a second 
order effect significant only during transient maneuvers when 
tracking is extremely difficult, this term is generally 
neglected. For the remainder of this report, the terms 
containing z and 6 will be dropped as in earlier LCOS mech¬ 
anizations. A later investigation is planned to determine 
the consequences of neglecting these terms. 

Equation (43) can be integrated on an on-board digital 
computer to obtain the lead ingle A(t). The speed of inte¬ 
gration can be modified to obtain better sight stability by 
time scaling equation (43). Let the new time variable t be 
defined such that 

t = Kt, dt = Kdx . (57) 

Substituting Eq (57) into Eq (43) and rearranging, we obtain 

T 



dA = K [W] )[Q]w 

eft 


V f A 


(D+T f D) 2(D+T f D) 


f [F x a N - (a N . A)F 


+ (a N . F)A]+ 


J v V a a. 


(D+DT f ) 


(58) 
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Equation (58) is a first order differential equation linear 

in X. The terms not containing X can be regarded as forcing 

terms, since they are dependent upon the attacker’s maneuvers. 

If steady state conditions exist, dX will approach zero, 

cTr 

and a steady state value of X can be obtained by integrating 

equation (44) until dX ~ 0, or, alternatively, dX can be set to 

cTr cTt 

zero, and X for the steady state situation can be obtained from 

the resulting algebraic equation. If steady state conditions 

do not exist, however, X(t) obtained by integrating equation 
(58) will lag the actual lead angle because of the time 

constant associated with X. 

This time constant can be modified merely by changing 
the time scaling factor K defined in equation (57). Alter¬ 
natively, K can be considered as a gain in a control system 
with output X(t). Increasing K decreased the time constant 
of X in equation (58); however, the greater K is, the greater 
will be the reticle sensitivity to changes in attacker 
aircraft acceleration, angle of attack, and rate of turn. 

Thus, a trade-off on the value of K is required between the 
pilot's ability to track a target with the reticle and the 
response of the reticle to changes in aircraft dynamics 
commanded by the pilot. The value K - .80 corresponds to 
the sight damping factor o = .25 (ref 1) found in most 
present day LCOS. 
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Transient Response Versus K 
The transient lesponse of equation (58) to a step 
input is illustrated in Figure 5 for various values of K. 
Initially, A(t ) = 0, and the step input of u(t) =* 0, 

i O 

_ A _ 

a^(t) = -*32.2k, a(t) * 0, t>0 corresponds to straight 
and level flight. Figure 5 clearly shows the lag existing 
in the LCOS in any nonsteady situation. The obvious 
result is to choose the gain K as high as possible 
consistent with the pilot's ability to track his target 
with the more responsive reticle. 


o 



TIME (SEC) 

Figure 5. Transient Response of Lead Angle for Various 
Values of K. 
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Sensitivity of X(t) to Accelerations in x and y 
The present LCOS measures only acceleration normal 
to the aircraft body x-y plane. In this report no assump¬ 
tion was made to measure acceleration only in the body z direc¬ 
tion; hence*it is of interest to determine the sensitivity of 

the lead angle Aft) to accelerations in the body x and y 
directions. A(t) for various accelerations*including the 

body x and y accelerations, is compared in Table 1 to the 
lead angle generated by assuming only normal acceleration. 

The results show that A(t) is quite insensitive to x while 
it is considerably more sensitive to y. In judging the 
importance of the sensitivity of A(t) to y,one must keep 
in mind that for coordinated flight y = 0, and for tracking 
situations it is most likely that the pilot will be coor¬ 
dinated. Therefore*LCOS can employ only normal acceleration 
z with little degradation in sight accuracy provided the 
pilot makes coordinated turns during tracking. 

Table 1. Sensitivity Of Lead Angle A(t) To Body x and y 
Acceleration 

Range = 2000' v = 800 fps h = 10000' 

Level Flight 


Acceleration 

2 


X (ft/sec ) 

y (ft/sec ) 

El 

16.1 

0 . 

0 . 

32.2 

0 . 

0 . 

0 . 

8.05 

1.2 

0 . 

16.1 

2.5 


Effect (Mils) 
Az 


.02 

.03 

0 . 

0 . 
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Further Simplification of Lead Angle Equations 
As stated earlier, the sensor information available and 
the geometry of the gun placement have an effect on the com¬ 
plexity of the terms in Eq (58). The sensors under discus¬ 
sion are angle-of-attack meter and body mounted accelerometers, 

Angle of Attack 

Generally the angle of attack of the aircraft* is 
measured in the body x-z plane; hence*for small angles it 
can be considered as a rotation about the body y axis. The 
angle of attack of the aircraft can then be expressed as 

(assuming no side slip) 


5= (:)- 


(59) 


The angle of attack of the gun <* G is the angle required in 
Eq (58), and the relationship between o and a G is 

a G = a - F g (60) 

where Fg is the angle considered as a vector that the gun 
makes with the x axis of the aircraft. Figure 6 shows the 
geometrical relationship between a G , a, and Fg. If the 
gun is mounted in the body x-z plane, then 



x^ (aircraft x axis) 

F (gunline) 

» 

va (velocity vector) 


Figure 6. Geometrical relationship between a r , a, and 0. 


27 





















a G 15 
a- - j. ; 


Accelerometers 

If only one accelerometer is employed in LCOS to measure 
body acceleration in the 2 axis, then a M in Eq (58) is 


r «' III = ji! 


and the terms containing simplify to 
F x a„ > -lb, 2 I 


a N - F = zb 3 


a N ' ^ ~ *X 3 • 


Gun Location 

If the gun is located in the body x-z plane, then the 
gun vector F becomes 


F = o 


I b, 

\ ' I 
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V. NUMERICAL COMPARISON WITH OTHER GUNSIGHTS 
A brief description of four other gunsights with which 
a comparison of AFALCOS will ^e made follows. A Fortran 
language listing of AFALCOS is found in Appendix A. 

•'HOTLINE" Historical Tracer 
This gunsight is a historical tracer type sight in 
which the relative trajectory of the bullets with respect to 
the attacking.aircraft is computed and displayed on a heads 
up display (HUD). If range is known from radar, a small 
reticle appears on the HUD at the point where a bullet at 
target range would be if it had been fired one time of flight 
ago - thus the name "historical tracer.” Very few assumptions 
are made in the calculations to obtain the relative bullet 
stream position; therefore, this sight is considered an 
accurate standard for comparison purposes. The authors are 
grateful to Honeywell Inc. for making this sight digital 
program available to the Air Force. The equations and pro¬ 
gram are considered proprietary by Honeywell Inc.; hence, 
they are not included in this report. 

Advanced Fighter LCOS 

A digital model representative of the lead computing 
sight for a current advanced fighter aircraft was programmed 
by members of the Astronautics and Computer Science Department 
of the US Air Force Academy. This digital model includes a 
digital representation of the lead computing gyroscope 
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employed in this lead computing sight. The digital gyroscope 
model was programmed by the Air Force Avionics Laboratory, 

WPAFB, Ohio. A listing of the Fortran language digital pro¬ 
gram is found in Appendix B. 

Honeywell Digital LCOS 

The equations and digital program for the above advanced 
fighter LCOS were given to Honeywell Inc. for development of 
a digital LCOS to be used in the Air Force Avionics Labora¬ 
tory's Comparative Gunsight Test Program in which a flight 
test comparison between a historical tracer type sight 
(HOTLINE) and an LCOS sight was to be made. Honeywell simpli¬ 
fied the equations and digital code to develop a compact code 
which could be programmed on an aircraft digital computer. 

It should be noted that during the course of this simplifica¬ 
tion and program checkout, Honeywell's engineers discovered 
the roll coupling discrepancy between the advanced fighter 
LCOS and their own HOTLINE and brought the matter to the 
attention of the Air Force. A Fortran language program list¬ 
ing is found in Appendix C. 

Air Force Avionics Laboratory LCOS 

This digital LCOS was developed by the Air Force Avionics 
Laboratory to support the laboratory's Digital Avionics Infor¬ 
mation System (DAIS) project. This development is similar to 
the other digital LCOS derivations except that the small angle 
assumptions made in the other LCOS derivations have been avoided. 
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Because this LCOS uses a different time of flight equation than 
the other LCOS's, to obtain a better basis for comparison the 
bullet time of flight calculation, was removed and the time of 
flight computed in the other LCOS’s was used. It should be 
noted that when the original time of flight calculation is 
used, worse agreement is obtained between this LCOS and HOT¬ 
LINE. A Fortran language listing is found in Appendix D. The 
development of the equations for this LCOS is found in 
reference [2]. 

Comparison 

Five cases at 2000 feet of range and five cases at 3000 
feet of range were used as a basis of comparison. In each 
case the altitude was 10000 feet, aircraft angle of attack 
a was zero, gun angle 0 G was zero, muzzle velocity was 3300 
ft/sec, and sight damping factor a for the Honeywell, advanced 
fighter, and AFAL LCOS's was 0.25. Aircraft load factor in 
each case was 4. The sights were allowed to settle for 12 
seconds so that a steady state comparison could be made. 

Table 2 shows the results of the comparison at 2000* range, 
whereas Table 3 contains the comparisons at 3000' range. 

Table 4 contains the results of the comparison between HOTLINE 
and the AFAL LCOS at both 2000' and 3000' of range. The time 
of flight calculation used in the AFAL LCOS is identical to 
that used in the other LCOS's. 
































































































































TABLE 4. Air Force Avionics Laboratory LCOS Compari 
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AAz and AElare the differences in mils between the HOTLINE 
sight bullet at target range as seen in a heads up display 
and the LCOS reticle position in the HUD. 6 is the radial 
error which is the square root of the sum of the squares of 
AAz and AE1. The Azimuth and elevation of the bullet at 
target range from the HOTLINE sight are listed in the columns 
under HOTLINE. The geometry is shown in Figure 7. 
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From Tables 2 and 4 where the data corresponds to a 
range of 2000* it is clear that for the first four cases the 
Honeywell LCOS, the AFAL LCOS, and AFALCOS are excellent 
while the advanced fighter sight is still very useable since 
the 3.4 mil error is well within the M-61 gun bullet 
dispersion pattern. When a high roll rate (p = .2 rad/sec) 
is introduced, however, the Honeywell LCOS, AFAL LCOS and 
advanced fighter sight introduce large errors--17.7 mils, 

12.0 mils, and 21.5 mils respectively--because of the 
incorrect mechanization of the roll coupling in these sights. 
The AFALCOS for this case has an error of only 2.2 mils--an 
improvement by a factor of ten. 

When the range is extended to 3000’ for the first four 
cases, the LCOS sights degrade to an average radial error of 
6.2 mils for Honeywell LCOS, 9.2 mils for AFAL LCOS, 10.2 mils 
for the advanced fighter LCOS, and 5.2 mils for AFALCOS. 

When a roll rate of 0.2 rad/sec is introduced, the 
Honeywell LCOS, AFAL LCOS, and advanced fighter LCOS have an 
average radial error of 56.0 mils, 38.4 mils, and 66.8 mils 
respectively, while AFALCOS has a radial error of only 5.7 
mils--again, an improvement by a factor of ten. 
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VI. CONCLUSIONS 


(1) This report presents the derivation of the equations 
necessary to implement a lead computing optical sight in a 
fighter type aircraft- A fortran digital program is 
developed to implement these equations on a digital computer. 

(2) The equations developed require the angular rate 
vector of the aircraft including the roll rate. These signals 
may be obtained from either a three axis strapdown gyroscope 
system or an inertial navigation system employing a stable 
platform. Use of either of these angular rate measuring 
devices will permit the elimination of the two axis lead 
computing gyroscope presently employed in LCOS. 

(3) The assumptions leading to incorrect roll coupling 

in the Honeywell and the advanced fighter LCOS sights have 

% 

been eliminated, and the roll coupling in AFALCOS now corre¬ 
sponds very closely in the steady state with the results of 
the HOTLINE historical tracer sight. 

(4) In the comparisons made, in Tables 2, 3 and 4 the 
AFALCOS steady state radial error is less in almost every 
case than that of the Honeywell, advanced fighter and 
Avionics Lab LCOS sights. Therefore, when the HOTLINE sight 
is used as a steady state accuracy standard, AFALCOS is 
judged to be a more accurate sight than any other LCOS known 
to the authors. 










(5) The digital code is not extensive for AFALCOS; 
hence, no particular programming difficulties are antici¬ 
pated for airborne application. 

(6) A word of caution--the numerical comparisons made 

in section V are steady state comparisons. Before final 
suitability of any lead computing sight can be established, 

the entire system including pilot, target, sight display, 

flight control system, and aerodynamic response of the 

aircraft, must be tested in a dynamic environment. 
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APPENDIX A 


This appendix contains a listing of the Fortran 

digital code for AFALCOS. The ?nput/output is structured 
to operate on a remote terminal of the Burroughs B6700 

digital computer at the U. S. Air Force Academy. 
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C- LEATHAM,A.L. USER CODE BYZ0185. DFACS. EXT. 2136 

SRESET FREE 

FILE 5=INPUT2,UNITsDISK,BL0CKING=30,RECORDS 14 

FILE 7=LINE,UNIT=PRINTER 

FILE 6sOUTPUT,UNIT=REMOTE,RECORDs2 2 

1000 FORMAT(1H0,F4.2,3F10.3,2F10.2) 

1001 F0RMAT(1H0,'TIME , .6X •ELA* t 7X, , ALA , ,5X,*FLITlME', 

S2X,'ELHUO(MILS)*,1X, AZHUD{MILS)',//) 

REAL KB.JV 

NAMELIST/INPUT/P,Q,R,D,D00T,VA,AX,AY,AZ,AL,GA,GAIN,DT,FTIME,HA 
S ,VM,IPRINT 
C 

C INPUT DATA 
C 

C P,Q,R ... BODY ANGULAR RATES IN RADIANS/SEC 
C D,DD0T — RANGE AND RANGE RATE IN FT AND FT/SEC 
C VA -- AIRCRAFT SPEEO IN FT /SEC 

C AX,AY,AZ- SOOY ACCELEROMETER MEASURED ACCELERATION IN FT/SEC**2 
C AL.GA — AIRCRAFT ANGLE OF ATTACK AND GUN ANGLE IN RADIANS 
C GAIN SIGHT SENSITIVITY PARAMETER .8 NOMINALLY 

C DT.FTIME- INTEGRATION STEP SIZE AND FINAL TIME IN SEC. 

C KA,VM — - ALTITUDE IN FEET AND MUZZLE SPEEO IN FT/SEC. 

C IPRINT —- PRINT CONTROL INDEX 
C 

500 READ(5, INPUT) 

PRINT INPUT 
C MCD ATMOSPHERE 

IF(HA.GT.36000.) GO TO 50 

RHQ=(,034475+(.019213E-10*HA-.050331E-5)*HA)**2*2, 

Ssll 17. 1-.00412778*HA 
GO TO 55 

50 DHsHA-36000. 

RH0=(.018828+(.039227E-10*DH-.043877E-5)*0H)**2*2. 

S = 968 5 

55 ORATIQsRHO/.00238 


a 





s 

K 

i 

11 


VPsVM+VA 
KBs 00614 

VLS=D*KB*VP*<\5*DRAT10 

VCs-DOOT 

VOSsVM+VC-VLS 

VCM=(VOS**2-4.*(VA-VC)*VLS)**.5 

RTF=,5*(V0S*YCM)/D 
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; TF IS 8ULLET TIME OF FLIGHT 
TF=t./RTF 
VF=9/TF«VC 
JV=(VM-VF)/(VA+VM) 

; 8 IS THE GUNLINE UNIT VECTOR IN BODY COORDINATES 

B1=COS(GA) 

B2=0. 

B3=SIN(GA) 

AL IS THE A/C ANGLE OF ATTACK IN RAOIANS* PCS. FOR POS. LOAD FACTOR 
GA IS THE GUN ANGLE MEASURED IN RADIANS DOWN FROM X AXIS 
GAL IS THE GUN ANGLE OF ATTACK 
GAL=AL-GA 
C6=D+0DOT*TF 
C7=TF*0*P/2./C6 
C1=VF/C6 
C2=JV*VA/C6 
C3=TF/2./C6 
BXAN2=-B!*AZ+AX*B3 
BXAN3=B1*AY-AX*B2 
PRINT 1001 

PRINT 1000,T.ELA,ALA,TF,ELAH,ALAH 
I SL IS THE SIGHT LINE UNIT VECTOR 
100 SL1=B1 + B2*ALA-B3*ELA 
SL2sB2-B1*ALA 
SL3=B3+81*ELA 

W2=-C1*ELA-C2*GAL*C 3*(BXAN2-< AY*ELA +AZ*ALA)*B2+(AX*B1♦AY*B2+AZ* B3) 
3 *ELA) 

W3=-C1*ALA+C3*(8XAN3-(AY*ELA+AZ*ALA)*B3*(AX*BH-AY*B2+AZ*B3)*ALA) 
W2=W2+(-SL1*SL2+C7*SL3)*P+(1.-SL2**2)*Q-(SL2*SL3+C7*SL1)*R 
W3-W3-( SL1*SL3+C7*SL2)*P+(-SL2*SL3'*’C7 l5r SLI )*Q+( 1 .-SL3**2)*R 
C5=C7*SL1 

DETsI,-SL2**2-SL3**2^C5**2 

AINVW2=( (1 .-SL3**2)*W2*(SL2*SL3«-C5)*W3 )/DET 

AINVW3=((SL2*SL3-C5)*W2 + (1,-SL2**2)*W3)/DET 

DELA=GAIN*AINVW2 

DALA=GAIN*AINVW3 

; ELA AND ALA ARE EL. ANO AZ. LEAD ANGLE COMPONENTS W.R.T, GUN 

ela=ela*qt*oela 

ALA=ALA+DT *DALA 
T=T+DT 

ITIME = ITIME*H 

IF(ITIMEiEQ.IPRINT) PRINT 1000,T,ELA,A' *,TF,ELAH,ALAH 
IFUTIME.EQ.IPRINT) ITIME=0 

ELAB AND AlAB ARE EL AND AZ LEAD ANGLES IN BODY COORDINATES 
ELAB=B1*ELA 
ALAB=B1*ALA 

; ELAH ANO AZAH ARE HUD EL AND AZ ANGLES IN MILS. HUO ANGLE=0. 
ELAH--ULA-HUDZ/(VF*TF)+GA}*1000. 

ALAH=-(ALA-HUDY/(VF*TF))*1000. 

IF(T.LT.FTIME) GO TO 100 
ELAsO. 

ALAsO. 

T=0. 

GO TO 500 
END 
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APPENDIX B 

This appendix contains a listing of the Fortran digital 
code which represents the lead computing optical sight in 
an advanced fighter aircraft. The listing is a subroutine 
designed to operate in a larger program on the Burroughs 
B6700 computer at the U.S. Air Force Academy. 

SUBROUTINE F15LC0 

COHMON/HLLC OS/XHL.YHL,RHL.TQFHL 

COMMON/GUNOAT/H(I 50),S(300),TT(44),T 

COMMON/ANGLE/BPS I,BTHE,LEH,LAH,DELXX,OELYY,ELAH, ALAH.GAIN 
COMMON/SGYRQ/AG,EG.SEG,CEG.SAG,CAG,DT 

COMMON/ICS/OX.DY,DZ,BKORHO,BO,TAU,R1500,DELTA ,PCY,Ptt, 

.ALCDFI.THE,PSI. ALCDF.DELT.THG.ALOQP.RHO»VM,AT F 
REAL LAH.LAG.LEG.LAC.LEC.LAH.LEH 
EQUIVALENCE 

.(HA.S(24)),(VTAS,S(3t)),(R .S(8)), 

,(VC,S(9)),(SAL,S(35)).(SBE ,S(36)), 

.(AXB.S(5)),(AYB,S(6)),(SA ,S(7)), 

. (PP,S(2)),(QQ,S(3)),(RR,$(4)), 

.(ALPHA.S(201>),(BE.S( 202)) 

NAMELIST/DATTA/VTAS.R,HA,ALPHA,SA.SJ.PP.QQ.RR 
.,XRI,YRI.ZRI,XDOTP,YDOTP.HDOTP,XDOTE,YDOTE,HOOTE 
NAPE LI ST/OUTPUT/RHO,VLS,VOS,VCH.RTF,XI,VR,VCK,P.LAH,FG. 
.LEG.LAG.XH.YH.ZH, CEG,CAG,SEG.SAG 




i 

IF(T.GT.0.06) GO To 89 

s. 

i 

S' 

k 

C IMTRALIZE SINES & COSINES 

| 

i 

CAG=1. 

| 

e 

s 

SAG=0. 


1 

1 

SEG=0. 



CEG=1. 


1 

I' 

AG=0. 

& 

EG=0. 


1 

891 CONTINUE 

£ 

V 

l 

SSA=-SA 

1 

i 

ALPHH=ALPHA 

1 

t 

1 

RHD=RHO/.00238 

S' 

i 

4. 

BK=0,00614 

%' 

l 

I 

C VELOCITY OF HISSILE(FT/SEC 

Ss 

VM=3 300 


1 

TAU=.25 

| 

* 

CGUN ANGLE 

*r 

? 

GA=0,0 

V 

U 

£ 

i 

i 

TGHE=0.0 

& 

% 
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C TGHH IS THE SUM OF GUN ANGLE AND HUD ANGLE. 
TGHH=TGHE+GA 
SK=.714236 
0L=0.0 
DA=0.0 
DE=0.0 
TGHA=0.0 

C COMPUTER CYCLE TIME IS DT 
DTs.06 

223 CONTINUE 

IF(VTAS ,GT, 2000 .OR. VTAS .LE. 0) VTAS=970 
C COMPUTE INITIAL VELOCITY IN AIR MASS 
1 VP=VM+VTAS 

C AVERAGE VELOCITY LOST FOR STATIONARY TARGET. 
VLS=R*BK*SQRT(VP)*RHD 

C AVERAGE OVERTAKING VELOCTIY FOR STATIONARY TARGETS. 
VC=0. 

VOS=VM*VC-VLS 

C VELOCITY CORRECTION FOR MOVING TARGET 
VCM=SQRT(VOS**2-4*<VTAS-VC)*VLS) 

C RECIPROCAL OF TIME OF FLIGHT CALCULATION 
RTF=0.5*(V0S»VCM)/R 

C RECRIPROCAL OF SIGHT SENSITIVITY. 

RTN=RTF>TAU*VC*(RTF+VLS/R)/VCM 
C OOME MAG CURRENT 
XI=RTN*SK 

C PROJECTILE AVERAGE RELATIVE VELOCITY. 

VR=RTF*R-VC 

C BALLISTIC CURVATURE coefficient. 

VCK-VTAS*<1-VM/VR)/VP 
C TOTAL GUN ANGLE OF ATTACK. 

ALPHAG=GAfALPHH 
C PRECESSION RATE COMMAND. 

P=.5*(SSA+R*SJ/VTAS)/VR«-RTN*VCK*ALPHAG 
222 CONTINUE 
100 CALL GSM(XI.P) 

C TOTAL LEAD ANGLE. 

LAM=ARCOS(CEG*CAG)*SK 
IF(ABS(SAG) .LE.1E-4) GO TO 110 
C LEAD ANGLE ROTATION. 

FG=ATAN2(SEG/CEG» SAG) 

GO TO 120 
110 FG=3.1415926/2 
C GYRO AZIMUTH LEAD ANGLE. 

120 LAG=ATAN<COS(FG)*TAN(LAM)) 

IF (ABS(LAG), LE .1 ,t-4) GO TO 130 
C GYRO ELEVATION LEAD ANGLE, 

LEG=ATAN(TAN<FG)*SIN(LAG)) 

GO TO 140 
130 LEG=LAM 

IF (EG.LT.O.) LEG=-LAM 
C CORRECTED ELEVATION LEAO ANGLE. 

140 LEC=ATAN((1*DL/R)*TAN(LEG) -DE/R) 
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♦ DA/R) ♦TGHA 


CORRECTED AZIMUTH LEAO ANGLE. 

LAC=ATAN((1♦Dl/R)*TAN(LAG) 

HUD COORDINATES OF LEAD ANGLE 

XH=C0S(LEC)*C0S(LAC)*C0S(TGHE)+5IN(LEC)*SIN(TGHE) 
YH=COS(LEC)*SIN(LAC) 

ZH=-COS(LEC)*COS(LAC)*SIN(TGHH)♦$IN< LEC )*COS<TGHH) 
HUD AZIMUTH LEAD ANGLE COMPONENT. 

LAH =1000.*ATAN (YH/XH) 

HUD ELEVATION LEAO ANGLE COMPONENT. 

LEH =1000.*ATAN (ZH/XH) 

998 CONTINUE 

999 CONTINUE 
0ELXX=LAH-8PSI 
DELYY=LEH-BTHE 
CONTINUE 
RETURN 

END 

C LCG GYRO SIMULATION MODEL 
C 

CINPUTS 


303 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 


MAG CURRENT COMMAND 
PRECESSION RATE COMMAND 
ACFT ROLL RATE 
ACFT PITCH RATE 
ACFT YAW RATE 
OUTPUTS 

AG - GYRO AZIMUTH GIMBAL ANGLE 
EG GYRO ELEVATION GIMBAL ANGLE 
DEFINITIONS 
SG -GYRO SPACE ANGLE 


XI 

XP 

P 

Q 

R 


TG -GYRO TIME CONSTANT(SENSIVITY) 


TM - MOO IF I ED BYRO TIME CONSTANT 
TGL- GUNLINE ELEVATION ANGLE 
XG - GYRO ROLL RATE 
YG - GYRO PITCH RATE 
ZG - GYRO YAW RATE 
OT - COMPUTER CYCLE TIME 
SUBROUTINE GSM(XI.SP) 

CCMMON/GUNDAT/H(150).S(300),TT(44),T 

COMMON/ANGLE/BPS I,BTHE,LEH,LAH,DELXX,DELYY,ALAH,ELAH,GAIN 
CQMMON/SGYRQ/AG,EG,SEG,CEG,SAG,CAG,DT 

COMMON/IC S/DX,DY,DZ,BKORHO, BO ,T AU, R1500,DELTA ,PCY,PCZ, 

.ALCDFI ,THE,PSI, ALCOF,DELT,THG,ALOOP,RHO,VM,ATF 
REAL LAM,LAG,LEG,LAC,LEC,LAH,LEH 
EQUIVALENCE • 

,(HA,$(24)),(VTAS,S(31))»(R ,S(8)), 

,(VC,S(9)),(SAL,S(35)),(SBE ,S(36)), 

,(AXB,S(5)),(AYB,S(6)),(SA ,S(7)), 

. (PP,S(2)),(QQ,S(3)),(RR,S(4)), 

.(ALPHA,S(201)),(BE,S(202)) 

P Q R IN RAD/SEC 

* *NAMELiST/OUTT/XG,YG,ZG,SG,TG,SM r EGD,AGD,EG,AG,SAG,CAG,SEG,CEG 

TGL=0. 

XG=PP*COS(TGL)-RR*SIN(TGL) 
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YG=QQ 

ZG=PP*SIN(TGL)♦RR*COS(TGL) 

SG=ARCOS(CAG*CEG) 

IF (A8S(SG),LT,1,E-4) GO TO 80 
TG=(SG/1.4)/(XI*SlN(SG/1.4)) 

SM=SG/SIN(SG) 

90 EGD=XG*SIN(AG)-YG*C0S(AG)-(C0S(AG)*$IN(EG)/1.0)*<SM/TG)-SP 
AGD=-ZG-TAN(EG)*(XG*COS(AG)« , YG v SIN(£G)) - 
. <SIN(AG)/(COS(EG)*SIN(SG)))*(SG/TG) 

GO TO 100 

80 EGD=XG*SIN(AG)-YG*COS(AG)-SP 

AGO =- ZG -T AN (EG ) * (XG*C OS (AG) ♦ ' Y G*S I N ( AG )) 

100 EG=EG+EG0*0T 
AG=AG*AGQ*DT 
$EG=SEG*EGD*OT*CEG 
CEG=CEG-EGD*OT*SEG 
SAG=SAG*AGD*DT*CAG 
CAG=CAG-AGO*DT*S AG 
RETURN 
END 
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APPENDIX C 

This appendix contains a listing of the Fortran 
digital code for the Honeywell digital LCOS. The 
listing is a subroutine designed to operate in a larger 
program on the Burroughs B6700 computer at the U.S. Air 
Force Academy. 

SUBROUTINE LCOSS 

cqmmon/hllcos/xhl,yhl,rhl,tofhl 

COMMON/GUNDAT/H(150),S(300)*TT(44),T 

COMMON/ANGLE/BPS I,BTHE,LEH,LAH,DELXX.OELYY,ElAH,ALAH.GAIN 
COMMON/SGYRO/AG,EG,SEG,CEG.SAG,CAG,DT 

COMMON/ICS/OX,OY,DZ,BKDRHO,BO,TAU,R1500,DELTA ,PCY,PCZ, 

.AlCDFI ,THE,PS!, alcdf,oelt,thg,aloqp,rho,vm,atf 

REAL LAM.LAG,LEG,LAC,LEC,LAH,LEH 
EQUIVALENCE 

.(HA,S(24)),(V ,5(31)),(RPA,S(8)), 

.(VC,S(9)),(SAL,S(35)),(SBE ,S(36)), 

,(AXB,S(5)),(AY8,S(6)), (ACGZ, S( 7)), 

. (P ,S(2)),(Q ,S(3)),(R ,S(4)), 

.(AL,S(20I)),(BE,S(202)) 

IF (T.GT..06) GO TO 500 

THG IS THE GUN ANGLE OF ATTACK 

THG=0. 

OX,0Y,OZ ARE DISTANCES FROM GUN TO HUD 
DX=0. 

DY=0, 

DZ=0, 

C BK IS AN EMPERICALLY DERIVED ALLISTJC COEFFICIENT. 

C BKDRHO : BK/RHO = .00614/.00238, 

BKDRHO=2,58 
BO:,003713 

C TAU IS THE AMPING FACTOR, 

TAU:.4 

TAU:.25 

R1500 : 1500. 

C R1500 IS THE FIXED RANGE WITH NO RADAR LOCK. 

C DELTA IS THE COMPUTER CYCLE TIME. 

PCY =0.0 
PCZsO.O 
THE : 0. 

PSI : 0. 

BPS 1:0, 

211 BTHE-O. 

ALOOP=50, 

0ELT:3. 

200 OELTAsDELT/ALOOP 
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P yz^m-v* 
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500 


VM IS THE MUZZLE VELOCITY 
VM=3300. 

ALCDF=0«*TAU)O,4 
ALCDF=1.25 

ALCOFI s ( J/O+TAU) ) 

alcofui./alcof ; ; 
continue 

Af||=0 INDICATES RADAR LOCK 
A393=0, 

XJJf I s T «£ TRUE AIR SPEEO. 
VTAS=V 

RANGE=R?500 

JT<A393 .EQ.O.) RANGEsRPA 

Cf _ iE5.5f A * LT • J •> RANGErl. 

510 CONTINUE 




RHO VS RANGE FOR CHECK. 


THIS IS A VALIDITY CHECK, nnu v 

B=BO^'.000000888*(1 255.-(V-VC))) 

SAVE s-.000000618*RANGE♦6 
IF(RHQ.LE.SAVE) GO TO 556 
RTN-,33333 
GO TO 557 

556 CONTINUE 
VP=VM*VTA$ 

VP IS INITAL PROJECTILE VELOCITY IN AIRMANS 
VLS=RANGE*8KDRH0*SQRT(VP)*RHO 
VC =0, 

VOS=Vm+VC-VLS 

CHECK=V0S**2-4,*(VTAS-VC)*VLS 
VCM=SQRT(CHECK) 

VF=.5*(V0S+VCM) 

mfM !* 8 PROJ£ct,u VEL0C,Ty - 

RT ^i S RECIPROCAL OF time of flight. 

atf=i,/rtf 
VR=RTF*RANGE-VC 
TFDOT=VC*<RTF*VLS/RANGE)/VCH 
TEOOT IS the RATE OF CHANGE OF T0F/TF. 

rtn=rtf«tau*tfoot 

557 TRTni|jRTN RECRIPR0CAL 0F SIGHT SENSIIIV,TY - 
TRTN IS THE SIGHT SENSITIVITY, 

TF = TRTN 
VT s VM ♦ V 

VVT IS T rv/VT TAL IN,TI4L WOJECTIIE VELOCITY IN AIRMASS 

ALCPSI=ALCOF*PSI 
QXYsQ-P*ALCPSi 

RZG=R+ALCDF*THE*(P+Q*ALCPSI) 

RI=(VM-VR)/VR 

ALG =RI*VVT 

ALGI = (AL+THG)*ALG 

TFI2=.5*TF 
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PCY=RZG*TF 
ACCZ= ACGZ 

PCZ=QXY*TF -ACCZ*TFI2/VR-ALG1 
XCAl=ALCDF1*0E LTA/T F 
THE=THE-< THE+PCZ)*XCAL 
PSI=PSI-(PSI«-PCY)*XCAL 
BPS! s PS! ♦ OY/RANGE 

BTHE = THE - DZ/RANGE 

*-THG 

BPS!=1OOO.*BPS I 

BTHE=!000.*BTHE 

RETURN 

END 






APPENDIX D 


This appendix contains a listing of the Fortran digital 
code for the Air Force Avionics Laboratory digital LCOS. 

The code was extracted from reference [2]. The program is 
designed to be run from a remote terminal of the Burroughs 
B6700 computer at the U.S. Air Force Academy. 


$RESET FREE 

FILE 5=AVDATA,UNIT=DISK,BL0CKING=30,REC0RD=14 
FILE 6=QUTPUT,UNIT=REM0TE,REC0RD=22 

REAL KBRHO,KH,KSIG,LE,LA, LAD, LED 
DIMENSION TME<205),PLE(205),PLA(205),IBUFF(1024) 

1 IPT=0 

C THE FOLLOWING DATA REPRESENT INPUT VARIABLES 
C CONSTANT INPUTS ARE ASSUMED HERE AS AN EXAMPLE 

namelist/qata/range,vc.p,q,r,an,alpha,va,rhq,ha 
REA0(S(DATA) 

PRINT DATA 
C RANGE IS IN FT 

C VC IS NEGATIVE OF RANGE RATE IN FT/SEC 

C P,Q,R ARE ANGULAR RATES OF AIRCRAFT IN BODY AXES, RAD/SEC 
C AN IS NORMAL ACCELERATION IN FT/SEC/SEC (-32.17 WHEN STRAIGHT & LEVEL) 
C ALPHA IS ANGLE OF ATTACK IN RADIANS 
C VA IS TRUE AIRSPEED IN FT/SEC 
C RHO IS AIR DENSITY IN SLUGS PER CUBIC FOOT 
TIME=0. 

C SET CONSTANTS 
C COMPUTER CYCLE TIME IS DT 
DT=.06 

C SIG IS THE SIGHT OAMPING FACTOR 
SIG=0.4 

KSIG=1./(1.*S!G) 

C VM IS MUZZLE VELOCITY IN FT/SEC 
VM=3300. 

C KBRHO IS A BALLISTIC COEFFICIENT DIVIDED BY SEA LEVEL AIR DENSITY 
KBRH0=0.00614/.00238 

C GA IS GUN ANGLE WITH RESPECT TO X AXIS OF AIRCRAFT 
GA=0.0 
CQSGA=1.0 
$INGA=0.0 

C RRH IS RECIPROCAL OF GUN HARMONIZATION RANGE 
RRH=.0005 

C DA AND DE ARE Y AND Z DISTANCES OF GUN FROM HUD IN FT 
DA=0. 

DE=0. 

CINITIALIZE 
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C IPRINT IS A PRINT CONTROL INDEX SET TO PRINT EVERY 5TH TIME INTERVAL 
IPRINT=0 
LA=C. 

LE=0. 

R0C=O. 

SQRV = S QRT (VA+V M) 

VLS=KBRHO*RHO*RANGE*SQRV 

VQS=VM+VC-VLS 

VCM=SQRT(VQS*V0S-4.*(VA-VC )*VLS> 

PRINT 10 

10 F0RMAT<1H0,6X,'TIME',6X, , LA\7X,»LES7X # , LAD\7X,*LED , ,7X, , WK\ 

37X,»WJ*//) 

C END INITIALIZATION — BEGIN LOOP 
DO 3 1=1*200 
IPT=IPT+1 

SLA=LA*(1.-.16667* LA* LA) 

SLE=LE*( 1. 16667* LE* LE) 

CLE=1 «-»5*LE* LE 
CLA=1 .-.5* LA*LA 
RC0S=1./(CLA*CLE) 

RDE=(-VC-RDC*SLE)*RC OS 
RE=RANGE*RCOS 

C APPROXIMATION USED FOR SQUARE ROOT 

C IF Y=APPROX. SQUARE ROOT OF X, THEN SQRT(X>=.5*<Y+X/Y) IS VERY CLOSE 
SQRV=.5*(SQRV+(VA*VM)/SQRV) 

VLS=KBRHO*RHO*RE*SQRV 

VQS=VM-RDE-VLS 

C APPROXIMATION USED FOR SQUARE ROOT 

C RDE=-VC MAKES TOF CALCULATION IDENTICAL TO OTHER LCOSS. 

RDE=-VC 

VCM=.5*(VCM^(V0S*V0S-4.*(VA+RDE)*VLS)/VCM) 

VE=«5*C VOS+VCM) 

C RTF ANO VF ARE 1/TIME OF FLIGHT AND AVG RELATIVE VELOCITY OF BULLET 
RTF=VE/RE 
VF=VE+ROE 
RRANGE=1./RANGE 
KH=1.-RANGE*RRH 

VN=VF-SIG*RDL VLS*(VF*VA)/(VCM*VF) 

PG=P*COSGA-R*SINGA 

RG=P*SINGA+R*COSGA 

RDC = VA*(ALPHA*GA)*(VM«VF)/(VA*VM)+« 5*AN/RTF 
C WJ AND WK ARE COMPUTED SIGHT LINE ANGULAR RATES 
WK=(KH*OA*CLA^RTF-VN*SLA)*RRANGE 
WJ=((ROC-KH*DE*RTF)*CLE-VN>CLA*SLE)*RRANGE 
LED=(WJ*SLA*PG-CLA*Q)*KS IG 
LAD=((WK-SLE*(CLA*PG«’SLA*Q))/CLE-RU>*KSIG 
LA=LA+LAD*DT 
LE=LE«-LED*DT 
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C LA AND LE ARE AZIMUTH AND ELEVATION ANGLES OF COMPUTED SIGHT LINE 
C LA AND LE ARE IN RADIANS. COMPUTE PLA AND PLE AS NEGATIVE IN MILS 
PLA(IPT)=-1000.*LA 
PLE(IPT)“-IOOO.*LE 
TIME=TtME+DT 
TME(IPT)=TIME 
IPRINTsIPRINT-M 
IF( IPRINT.NE.25) GO TO 3 
PRINT 11.TIME,LA,LE,LAU,LED,WK,WJ 
IPRINT=0 

It FQRMAT(1X,F 10.3»4F10.6,2F10.6) 

3 CONTINUE 
T0F=1./RTF 

NAMELIST/DEBUG/TOF, VF, VE, VCM, VOS, VLS 
PRINT DEBUG 
GO TO 1 
END 


SO 
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